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ABSTRACT 

We constrain the behavior of the radio luminosity function (RLF) of fwo classes of ac¬ 
tive galactic nuclei (AGN) namely AGN of low radio power (LRP) and BL Lac objecfs. 
The extrapolafion of the observed steep RTFs to low power predicts a space density 
of such objects that exceeds that of the sources that can harbor them and this re¬ 
quires a break to a shallower slope. For LRP AGN we obtain Pbr.LRP ^ 10^° ® W Hz“^ 
at 1.4 GHz to limit their density to be smaller than that of elliptical galaxies with 
black hole masses Mbh > By combining this value with the limit derived 

by the observations the break must occur at Pbr,LRP ~ 10^°'® — 10^^ ® WHz“^. For 
BL Lacs we find Pbr,BLLAC ^ WHz“^ otherwise they would outnumber the 

density of weak-lined and compact radio sources, while the observations indicate 
Fbr,BLLAC ^ WHz“^. In the framework of the AGN unified model a low lu¬ 

minosity break in the RLF of LRP AGN must correspond to a break in the RLF of 
BL Lacs. The ratio between Pbr,LRP and Pbr.BLLAC is 10^, as expected for a jef 
Doppler factor of ^ 10. 
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1 INTRODUCTION 

Many astrophysical quantities obey a distribution well de¬ 
scribed by power-laws. Important insights into the physical 
processes producing such distributions can be obtained from 
their power-law index. In addition, the location at which the 
distribution departs from a power-law is also of great im¬ 
portance. In most cases, a change in the distribution laws is 
observed at their high end, in the form of either a steepen¬ 
ing of the distribution, or an exponential cutoff. For exam¬ 
ple, the peak frequency in the synchrotron emission of active 
galaxies is associated to the largest energy attainable by rel¬ 
ativistic electrons and it carries essential constraints on the 
acceleration and cooling mechanisms. Equally important is 
the determination of the behavior of the distributions at the 
low end as it might reveal the presence of, e.g., a minimum 
mass for the formation of a star or a galaxy. Unfortunately, 
the low end behavior is often unaccessible to observations 
or, in other cases, it is blurred by the emergence of strong 
selection biases. 

In this Letter we focus on the radio luminosity function 
(RLF) of two classes of active galaxies at low redshift, z < 0.1. 
The first is formed by low radio power (LRP) AGN: within 
the considered redshift limit, most of the radio emitting AGN 
are indeed objects of low power, log(Pr/ W Hz“^) < 24. The 
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second consists of BL Lac objects. According to the unified 
model for AGN, these two RTFs are expected to be connected 
with each other, as LRP AGN should represent the parent 
population of BL Lac objects. In both cases the RLFs are ob- 
servationally well defined over a very broad range of radio 
power. The RLF of LRP AGN also presents a break at high lu¬ 
minosity. Conversely, at their low luminosity end, both RLFs 
are well described by a pure power law by current observa¬ 
tions, without clear signs of a change in their slopes. In this 
Letter we show that it is nonetheless possible to constrain the 
power at which a low luminosity break in the RLF must oc¬ 
cur. 

Throughout the paper we adopt iLo = 70kms“^ Mpc”*^, 
Um = 0.30, and Ua = 0.7 


2 THE RADIO AGN LUMINOSITY FUNCTION 

The number density of AGN cannot exceed that of 
their potential hosts: this argument has b een used by 
iMauch & Sadler! 120071) and IC^taneo & Best! l2009l) to con¬ 
strain the low luminosity behavior of their RLF. They 
found that a break must occur at a radio luminosity 
log(Pr/ WHz“^) ~ 19.5 and 19.2, respectively. We use 
the same rationale, but including the information that their 
hosts, as shown in more detail below, are associated almost 
exclusively to very large black holes. 

The RLF of nearby galaxies extending to the lowest 
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luminosity has been obtained by iMauch & Sadleil i2007ll . 
It is derived from a sample of 7824 radio sources from 
fhe 1.4 GHz NVSS survey cross-correlated with the second 
data release of t h e 6 degree Field Galaxy Survey (6dFGS, 
llones et ^l2005h . iMauch & Sadleii separated radio sources 
associated to AGN from sfar-forming galaxies relying on 
fheir optical spectra. Their resulting RTF covers the range 
log(Pr/ WHz“^) = 20.4 — 26.4 and shows a steepening at 
high luminosities (above log(Pr/ WHz“^) = 24.59). 

iBest & HeckmanI ll2012h obtained an independent esti¬ 
mate of fhe AGN RLE This is slighfly shallower (reach¬ 
ing a luminosity log(Pr/ WHz“'^) ~ 22) but, since their 
sample is drawn from the SDSS, they could take advan¬ 
tage of fhe optical data to better characterize their proper¬ 
ties. The main result is that their hosts are almost exclu¬ 
sively massive early-type galaxies with black hole masses 
(derived from the s tellar velocity dispersion adopting the 
iTrernaine et al.l20 ^1aw~) l arger than log {Mbh/Mq)^ 7-St (see 
also Baldi & CapettillloT^l . 

The SDSS survey provides us with a highly com¬ 
plete sample that can be used to estimate the num¬ 
ber density of galaxies wifh large black holes. To fhis 
purpose, we analyze the 818,333 galaxies (MPA-JHU 
sample hereafter) in the value-added spectroscopic cat¬ 
alog produced by the Max Planck Institute for Astro¬ 
physics, and the Johns Hopkins University, and avail- 
able at http://www.mpa-qarchinq.m pq.de/SDSS/ 
(Brinchmann et ai]l2004l; iTremoriti et alj|2004l) . We consider 
all galaxies from redshift a = 0.05 to z = 0.1 to ensure a high 
level of completeness (see below) within the largest possible 
volume, which results in 0.053 Gpc^. 

We then select elliptical galaxies in the MPA-JHU 
sample by setti ng a threshold to the concentration index 
Cr ^ 2.6 (e.g . ferateva et alJ l200ll; iKauffmarm et alj l2003l: 
iBell et alj2003h and to their stellar velocity dispersion adopt¬ 
ing, conservatively, o-gtar ^ 140kms“^, corresponding to 
log(MBH/ATQ) > 7.5, finding -^38,000 objects. The result- 
ing sample has a high level o f completeness. According to 
iMontero-Dorta & Pradal j2009l) the completeness of fhe SDSS 
decreases with decreasing apparent magnitude, starting at ~ 
95% at the SDSS spectroscopic limit of r = 17.77, and be¬ 
ing still higher than 80% at r = 13.25. The vast majority 
of fhe LRP AGN hosts have a magnitu de in the range -23.5 
< Mk < -26.5 jMauch & Sadleill2007ll . This translates into 
r = 13.3-16.3 at z ^ 0.05 ( and r = 14 8-17 8 at 2 = 0.1) 
having adopted r — K = 3.0 iChang et al.l2006h . This implies 
that most of such galaxies in the selected redshift range are 
included in the MPA-JHU catalog. 

We now estimate how the total number of predicted 
LRP AGN varies depending on the behavior of the RLF 
at low power. A low luminosity break must occur in or¬ 
der to avoid the d ivergence of fhe LRP AGN number. 
The IMauch & Sadler! results suggest that it is located below 
log(Pr/ WHz“^) ~ 21.6, where the limited number of ob¬ 
served objects makes the RLF shape uncertain. We assume 
that the LRP AGN number density follows, below a break 
luminosity Pbr.LRP, a power law with an index m. We inves¬ 
tigate the effects of varying m from 0 (a flaf number count 
distribution) to infinity (equivalent to a sharp cutoff below 
Pbr.LRp), see the left panel of Fig. [1] We integrate numeri¬ 
cally the RLF starting from log(Pr/ WHz“^) ~ 18.4, a fac¬ 
tor of 100 below the observational limit: the total number 


of predicted LRP AGN grows rapidly at decreasing Pbr.LRP 
(Fig. [T] right panel). For m = 0 it exceeds the total num¬ 
ber of available hosts when log(Pbr .LHp/ WHz~^) ~ 21.8, 
a value that is inconsistent with the Mauch & Sadletl results 
even assuming that all massive galaxies host a LRP AGN. 
For larger m (and already for m = 1) the tension with the 
observations is eased, as we obtain log(Pbr,LRp/ WHz“^) ~ 
20.4-20.7. However, this requires an occupation fraction / of 
100%, i.e., all massive galaxies must harbor a radio-source 
with at least Pr ^ Pbr.LRp. The bivariate RLFs, obtained by 
IMauch & Sadler! splitting the sample into bins of infrared lu¬ 
minosity (their Fig. 16), indicate instead values in the range 
/ ~ 10% — 60% depending on the host magnitude. This im¬ 
plies that the values of Pbr.LRP derived above should be con¬ 
sidered as stricf lower limits. 

We conclude that: 1) the RLF of LRP AGN must break 
to a shallower slope for radio luminosities smaller fhan 
log(Pr/ WHz“^) ~ 20.5 — 21.5, and 2) the number density 
below Pbr.LRp must significantly decrease, leading to a peak 
in its distribution. 


3 THE BL LAC OBJECTS LUMINOSITY EUNCTION 

In this section we consider the RLF of BL Lacs af 1.4 GHz. 
We adopt the same line of reasoning followed above. In this 
case, the constraint on the RLF is derived imposing the re¬ 
quirement that the total number of BL Lacs does not exceed 
the total number of appropriate (see below) radio sources in 
the same volume. 

iBest & HeckmanI i2012ll provides us with a sample of 
18,286 AGN obtained by combining the SDSS DR7 MPA-JHU 
sample with the NVSS and the FIRST surveys, down to a flux 
density level of 5 mjy in fhe NVSS. We consider fhe same vol¬ 
ume as in Secf. 2, i.e., 0.05 ^ 2 ^ 0.1. We select the sources 
with 1) a rest-frame equivalent width of all emission lines 
smaller fhan 5 A (fhe sfan dard threshold for a BL Lac classi¬ 
fication, IStickeTet^llg^, and 2) a flux of the central radio 
component in the FIRST images larger tha n 5 mJy. The latter 
requ irement (based on the measurement of iBest & Heckm^ 
I 2 OI 2 I 1 is needed to isolate potential BL Lac objects accord¬ 
ing to their radio core emission while excluding, e.g., dou¬ 
ble lobed sources lacking of a core. The fhreshold is set at the 
same radio flux limit that defines the selected sample and it is 
included self-consistently in the analysis that follows. We are 
left with 524 potential BL Lac objects. We further check their 
radio morphology: indeed BL Lacs are only associated with 
specific radio structures on the kpc scale, being dominated 
by a compact radio core and, in some cases, showing also 
a one-sided jet or a halo jAntonucci & Ulvestadl[l98a) . We 
found that 49% of the selected sources are instead well de¬ 
fined Fanaroff-Riley type I or II (FRI or FR II) radio galaxies 
with sizes larger than 10-20 kpc. The 267 remaining objects 
are not necessarily all BL Lacs and certainly include intrud¬ 
ers, for example, FR I with size smaller than ~10 kpc that, 
due to the limited resolution of the FIRST images, are not 
recognized as such. However, with only a handful of excep¬ 
tions, the selected objects are all massive early-type galaxies, 
conforming wifh the standard LRP AGN hosts. 

On the other hand, there is the possibility that BL Lacs 
are so nuclearly dominated to be confused wifh stars in the 
SDSS images (and for this reason they might not be included 
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Figure 1. Left panel: AGN number density as a function of the radio power Pr at 1.4 GHz. The results of iMauch & Sadl^ 1200/1) are shown 
as black dots and their fit as a solid blue line. The fit has been extrapolated below different values of the break luminosity, Pbr.LRP/ with 
power-laws of different slopes m (red for m = 0, green for m = 1, and blue for m = oo). Right panel: total number of predicted LRP AGN 
in the SDSS-DR7 area with 0.05 ^ 2 : ^ 0.1 as a function of the low luminosity break in the RLF, estimated for the three values of m. The solid 
horizontal line represents the total number of possible radio-galaxies hosts in the same volume; the other horizontal lines allow for smaller 
occupation fractions. The thick blue segments highlight the ranges of acceptable Pbr.LRP values. The vertical line marks the upper limit to 
f\r,LRP set by the observations. 


in the MPA-JHU galaxies catalog). However, this possibil¬ 
ity can be disc arded. First of all, in the Roma-BZCAT cata¬ 
log of BL Lacs llMassaro et alJl200^ there are only 8 objects 
with 0.05:^ 2 0.1 and covered by the DR7. All of them 

have been selected as SDSS spectroscopic targets. Further¬ 
more, the brightest of these objects has a radio luminosity 
of 1.3 X lO^'^ W Hz“^. We estimated the corresponding non 
thermal contribution to the optical (r band) magnitude for 
similarly bright objects by assuming, conservatively, a radio- 
to-optic al spectral index of 0.5 , typical of high energy peaked 
BL Lacs iFadovani et al.l2003h . We obtain r = 17.5 at z = 0.05 
(and r =19.0 at 2 = 0.1), ~ 1-A magnitudes fainter than LRF 
AGN hosts. Only much brighter BL Lacs might be confused 
with stars. However, very luminous sources have a low space 
density and indeed no such object, that would have been eas¬ 
ily discovered, is known within the volume considered. 

We conclude that our estimate of 267 objects is a robust 
upp er limit to the BL Lacs number in our volume. 

iFadovani et al.l l2007ll derived the RLF of BL Lacs at 5 
GHz based on the Deep X-ray Radio Blazar Survey. They ob¬ 
tained a RLF consistent with a single power-law with a slope 
of 2.12 ± 0.16 in the range log(P,./ WHz-^) - 23.8-26.8. 
By adopting this RLF we estimate the predicted number of 
BL Lacs in the selected volume when varying the luminos¬ 
ity, Pbr.BLLAC/ at which the RLF breaks, see Fig. [3 In our 
calculation we assume a null radio spectral inde^j to con¬ 
vert from 5 to 1.4 GHz and take into account that objects 
with radio flux density smaller than 5 mjy are not included 
in our sample (this effect causes the flattening of log N for 
log(Pbr,BLLAc/WHz“^) < 23 in Fig.0. We conclude that the 
number of BL Lacs predicted by extrapolating their RLF to 


^ This assumption is based on the average r adio spectral index of 
BL Lacs with 2 < 0.1, i.e. < Or >= 0.03 IStickel et alj|l99lh . A 
steeper radio spectrum would slightly increase the derived limits on 
Pbr.BLLAC- 



Figure 2. Total number of predicted BL Lacs in the range 0.05 2 ^ 

0.1 in the SDSS/DR7 area as a function of the low luminosity break 
of their RLF, Pbr.BLLAC/ ^rid of the slope m of the number density 
below the break. The solid horizontal line represents the total num¬ 
ber of radio galaxies (weak lined and core dominated, see text for 
details) in the same volume; the dashed horizontal line corresponds 
to a BL Lacs percentage of the selected radio-sources of 10%. 


low luminosities exceeds the total number of possible hosts 
for log(Pbr.BLLAc/ WHz“’^) < 23.3-23.5, depending on the 
RLF slop e below the bre ak. On the other hand, the RLF data 
points of iFadovani et al.l exclude log(Pbr,BLLAc/ WHz“’^) > 
24.5. 

We conclude that the RLF of BL Lacs must break 
to a shallower slope for radio luminosities smaller than 
log(Pr/ WHz-i) - 23.3-24.5. 
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4 DISCUSSION AND CONCLUSIONS 

The results presented indicate that the RTFs at 1.4 GHz of 
both LRP AGN and BL Lacs must decrease their slope at, or 
just above, the low luminosity limit at which they are cur¬ 
rently determined. This might appear as a contrived coinci¬ 
dence, while, conversely, it is the natural consequence of the 
behavior of the RLF at low power. In fact, the number of ob¬ 
jects below the break that can be detected in a flux limited 
survey depends on the radio power as , where m is 

the number density slope for < Phr- The sfrong depen¬ 
dence on Pr implies that these objects are missing altogether 
in the currently available surveys (or at most found in a small 
number) and this prevents the RLF determination below Pbr. 

For LRP AGN we find that the number density below 
the break must significantly decrease producing a peak in 
the distribution. It is tempting to associate the presence of 
this peak to a minimal level of emission that is invariably 
associated with a supermassive black hole located at the cen¬ 
ter of a massive galaxy, which unavoidably accretes gas from 
the interstellar medium. The peak might instead correspond 
to a minimum level of accretion at which a radio jet can be 
produced. A better definition of the RLF at low luminosity 
is essential in this context: in the first case, we expect to find 
that all potential hosts harbor a LRP AGN, albeit of very low 
luminosity while, in the second alternative, t he occupation 
fractio n would never reach 100%. Note that iKimball et al.l 
bOllh found a similar dramatic fall of the number density 
from the analysis of a sample of 179 QSOs (Mi <-23) from 
the SDSS in the redshift range 0.2 < z < 0.3, for radio lumi¬ 
nosities below log(Pr/ W Hz“^) ~ 21 (at 6 GHz). 

Unfortunately, we have only sparse information 
on the radio properties of massive galaxies below 
log(Pr/ WHz“^) ~ 20 . These can be obtained with 

deep targeted observations of nearby galaxies. In the Virgo 
clusters there are two (out of 11 ) giant elliptical galaxies 
lacking of any sign of ra dio emission down to a limit of 
log(Pr/ WHz“^) ~ 18.6 jCapetti et al.ll200^ . despite their 
large mass (M* ~ 1O^^ ®M 0 ) and large black hole mass 
(Mbh ~ 1O®-1O® ®M 0 ). The small number statistics prevent 
us to derive any firm conclusion based on these observa¬ 
tions, but they apparently favor the interpretation that not 
all massive galaxies are able to produce a radio jet. 

The two RTFs considered in this Letter are expected to 
be linked with each other. Indeed, according to the AGN 
unified scheme, BL Lacs are the beamed version of low lu- 
minosity AGN (see e.g. lUrrv & Padovanil Il995l: iTadhunt^ 
1 20081) , i.e., BL Lacs are objects in which a highly relativis¬ 
tic jet is seen at an angle close to our line of sight. This ac¬ 
counts for their extreme properties, such as strong flux vari¬ 
ability and apparent superluminal motion of radio compo¬ 
nents. The effects of relativis tic beaming on the RLF was an- 


alyzed in a series of papers 

Urrv & Shaferlll984IUrrv et alJ 

ll99ll:IUrrv & Padovanilll99ll. 

19951). If the RLF of the parent 


population of unbeamed objects is described by a broken 
power law, the break will produce a change of slope in the 
BL Lac RLF at ~ x Pbr.LRP, where 5max is the maxi¬ 
mum value of the jet Doppler factor distributiorfl The ratio 


^ The Doppler factor 5 is defined as <5 = [r(l — /3 cos 0)] where 


between Pbr.LRP and Pbr.BLLAC is a factor ~ 10 ®, which sug¬ 
gests a value of the Doppler factor (imax ~ 10, remarkably 
consistent with the typical results obt ained from the obser¬ 
vations and the mode l predictions (e.g. lGhisellini et ^Il993l: 
ISavolainen et alJboiOh . This analysis implicitly assumes that 
all LRP AGN produce relativistic jets, a requirement that is 
not necessarily met by the objects at the very faint end of the 
RLF. The consistency between the estimates of the Doppler 
factor we just obtained with the theoretical expectations ap¬ 
pears to confirm the overall validity of fhis assumption. 

How can we improve our knowledge of the RLF at low 
power? For LRP AGN the lower limit to the break power, 
log(Pbr,LRp/ WHz“^) ~ 20.5, is not far from what is cur¬ 
rently accessible to radio observations of large sample of 
galaxies such as the NVSS. A factor of ~ 10 improvement in 
the flux threshold is already accessible to the JVLA and it will 
be sufficient to measure directly the RLF of LRP AGN below 
the break value we predicted. While waiting for such a sur¬ 
vey it is possible to ta ke advantage of the stacking technique 
(e.g. lWhite et alj2007h to measure the median radio luminos¬ 
ity of the population of massive elliptical galaxies at sub-mjy 
levels. In case the number density presents a peak in its dis¬ 
tribution, as our results suggest, this median luminosity is 
predicted to be located close to the peak value. 
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